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ABSTRACT 


Element abundances of cosmic rays Li through Si with energy 
above 0.8 GeV/amu were measured on a balloon borne instrument 
containing a total absorption ionization spectrometer. Statistical 
techniques were used to analyze the five measurements of each 
particle to determine its charge and energy . The technique allows 
a determination of systematic errors to be made. Corrections for 
Landau fluctuations, spark chamber inefficiency, and background 
particles were included. Comparison with other published results 
is made. Differences in the shape of the spectrum determined from 
measurements of different workers indicate that the absolute inten- 
sity is still known to only +15% between 2 and 10 GV/c rigidity. 
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INTRODUCTION 

In order to study the composition and energy spectra of multiply 
charged cosmic ray nulcei we have flown an ionization spectrometer on 
balloons 1 * 3 . In Che principal balloon flight with this instrument we 
found a dramatic difference in the differential spectra of iron and 
carbon and oxygen nuclei 3 . A similar difference has been observed by 
Juliusson 4 . Because of the profound implications of these measurements 5 ’ 6 
for the nature of the cosmic ray sources and because of some quantita- 
tive differences of these measurements with those of other workers 4 ’' > 
we have undertaken a study of the systematic effects which might have 
influenced our earlier results ♦ As we 2 and Webber' in his review 
article noted previously, our fluxes of carbon and oxygen were somewhat 
lower than those reported by other workers. We speculated that this 
might be due to back-scattering from the massive spectrometer used to 
measure the energy. 

These effects were checked using instruments on a balloon flight 
whose primary objective was a measurement of the cosmic ray electrons 8 . 

A schematic cross section of the instrument is shown in Figure 1. The 
spectrometer was sufficiently deep to measure electron energies (18 
radiation lengths), but was not deep enough to measure nucleon energies 
with good accuracy. It was sufficiently deep for nuclear interactions 
to take place, and so will produce backseat tering effects if they are 
important. The charge detection section of the experiment included two 
Cerenkov detectors and two plastic scintillators. This four fold charge 
detection scheme was an improved version of that used previously. 
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The particle trajectories were measured using a wire grid spark 
chamber. On this flight an 8 deck chamber was used in contrast to the 
4 deck chamber previously flown. A new more accurate trajectory deter- 
mining algorithm was used 9 which is able to eliminate spurious sparks 
from delta rays of the heavy nuclef 0 . This algorithm also can find 
multiple trajectories to look at the effects of interacting particles, 
atmospheric showers, etc. 

We shall show in this paper that this more sophistocated detector 
allows us to conclude that the largest uncertainty in our previous 
results was due to uncertainties in the efficiency of the spark chamber, 
and not due to backscattering . We present here our results on the flux 
and charge composition between Li and Si at a geomagnetic cutoff of 
2.9 Gv/c. The results are compared with those of other workers 11, 12 
from this same latitude. Systematic differences are discussed. There 
are discrepancies not only in intensity but also in the spectral shape 
between 2 and 10 Gv/c. 

EXPERIMENT 

The apparatus has been described by Silverberg 8 and was launched 
on a balloon from Cape Girardeau, Missouri on October 1, 1972. It 
drifted less than 200 miles from the launch site before landing. The 
geomagnetic cutoff is 2.9 Gv/c. 

The elements in the charge detection section of this experiment 
were similar to those previously flown except that a previously used 
Csl mosaic scintillator was replaced by a Cerenkov detector. The two 
plastic scintillators were 50cm x 50cm Pilot B and measured the ionization 
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loss rate of the particles (Si under 1,94 gm/cm 2 of matter). Two Cerenkov 
detectors were the other charge measuring elements. One was a piece 
of Pilot 425 (acrylic plastic with a wave shifter) viewed edge on through 
adiabatic light pipes by four photomultiplier tubes. The other was a 
piece of fused silica in a white box viewed by four PMT's. Because the 
flight was designed to look at electrons, the phototubes for both these 
Cerenkov detectors were linear only to Z ^ 16, and so no results above 
this charge were obtained. 

It was not possible to completely reconcile the response of the 
two Cerenkov detectors as a function of p. The edges of the fused silica 
crystal were blocked by support brackets, and even though the surface 
was roughened by sand blasting there apparently was still considerable 
total internal reflection at large angles of incidence. In the crystal 
with wave shifter, apparently less than % of the collected light is 
shifted and isotropized. The effects of light collection are complicated 
as the cone angle of the light varies with velocity and the internal 
reflection of this light depends on the incident angle of the particle. 
The net result is that the two detectors had different responses as a 
function of velocity and zenith angle. The differences were at the ± 10% 
level and have been ignored in this analysis. However, they made it 
impossible to unravel energy spectra from the Cerenkov response. Similar 
effects could be important in experiments in which a single C detector 
is used to determine velocity 11,13 . 

The tungsten modules shown in Figure 1 are made of 1 radiation 
length slabs of tungsten and plastic scintillators and they constitute 
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the energy spectrometer. It was designed for studying electron cascade 
showers, but at low energy (for p <; 0.9) its response was proportional 
to energy for heavy nuclei. In this energy range the energy 
deposit is dominated by ionization loss. This was verified by scatter 
plots of the Cerenkov and the spectrometer responses. At high energies 
(p > 0.9) the spectrometer is so thin that the response depends strongly 
on the location of the first interaction. This effect could be drastically 
reduced by placing constraints on the location of the first interaction* 
because of the consequent reduction in statistical significance of the 
data this analysis has not been done. 

| DATA ANALYSIS 

The data for each event consisted of a spark chamber readout and 
five pulse height measurements. The pulse heights, of two scintillators, 
two Cerenkov counters, and the summed spectrometer, were corrected by 
empirical calibration curves to put them in units of minimum ionizing 
muons. Using the trajectory as determined by the spark chamber, the 
two scintillator and two Cerenkov signals were normalized to vertical 
incidence. Non uniformities as a function of position in these four 
detectors were removed by simple maps. Each detector was divided into 
5x5= 25 squares and the average pulse height of oxygen in each 
square was found. A simple function over area that described these 
average pulse heights at the center of each square was created for each 
detector and was used as the mapping correction divisor. The detectors 
were adequately uniform for this technique to work; the uncorrected 
oxygen pulse heights were far enough removed from the nitrogen and neon 
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pulse heights and flourine is suffeciently rare that a "beam" of oxygen 
could be used. The four charge detectors were then normalized to the 
same pulse height for relativistic oxygen nuclei. 

Spark chamber trajectories of all events were required to exit the 
bottom of the spectrometer. The energy deposited in the spectrometer 
was found by summing the pulse heights in the several, layers of scintil- 
lator and multiplying by 8.17 MeV/muon as determined in a calibration 
run at Brookhaven 13 . 

The spectrometer was 180 gm/cm 9 thick (0.75 proton mean free paths) 
and so particles of range fa 180 gm/cm 3 were stopped by ionization loss. 
This energy was 620 MeV/amu for Be, 1100 MeV/amu for 0 and 1600 MeV/amu 
for Si. At higher energies where the energy deposit was dominated by 
the particle interactions, the pions produced, and the subsequent elec- 
tromagnetic cascades, the fluctuations were very large and depended 
strongly on where the particle interacted. The spectrometer was thus 
too thin to make spectral measurements. However, the spectrometer did 
help differentiate between particles of low charge and high velocity 
and those of higher charge and lower velocity. In this case the 
scintillator response turned up (like a relativistic increase) at 
high energy. This turnup increased as Z increased and is related to 
the non-linearity of response of the scintillator and to the high 
energy delta rays which can leave the core of saturation and can 
penetrate down to the scintillator from the matter in the detectors 
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above. The net result is to make the separation of adjacent charges 
more difficult above p =- 0.8. The spectrometer, even though not accur- 
ate enough to produce energy spectra, helps to resolve this problem. 

After mapping, the five pulse height measurements are analyzed to 
find the charge and energy of each particle. Each event accepted was 
plotted as a point in a five dimensional space. Concentrations of 
points occured where the average pulse heights for the various elements 
occurred. As the average pulse heights for each element vary with 
energy, the concentrations of points lie about lines in the 5-dimensional 
space. 

This analysis represents an extension of the technique described 
previously 2 ’ 1 4 ' 1 5 ’ 16 for discrete responses. A mathematically rigorous 
treatment of the technique has been given 17 . 

In order to find the response lines, we made two dimensional his- 
tograms or scatter plots as shown in Figures 2 and 3. The two scintill- 
ator pulse heights were, on the average, equal and the Cerenkov pulse 
heights had small enough variations that we treated their averages as 
equal. Thus, we had only to determine the concentration lines in three 
dimensions. The (Sj + ^)/2 vs. ((^ + C^)/2 plot in Figure 2 and the 
spectrometer energy vs. (Cj + C^)/2 plot in Figure 3 were sufficient 
for our needs. The lines and the standard deviations (as functions of 
Z and energy) were determined along each axis for all reasonably abun- 
dant elements. It was assumed that only photoelectron statistics 
contributed to the Cerenkov standard deviation; i.e., that knock-on 
electron fluctuations were of no concern 18 . The Cerenkov standard 
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deviation at the highest energies was obtained from plots such as the 
plot for oxygen nuclei given in Figure 3. At high energy the width of 
the vertical tail is solely due to Cerenkov fluctuations. The Cerenkov 
standard deviation at any other energy was then scaled like the square 
root of the Cerenkov signal. The large uncertainty in energy at high 
energies can be seen in Figure 3. 

A computer program was written which finds the distance of the 
point for each event from the nearest line in terms of standard devia- 
tions 17 . Each charge line was defined by twelve 5-dimensional points 
and by linear interpolations between them. Interpolations were also 
made between charge lines to produce non- integral charge lines. (Charge 
was allowed to vary continuously and a response surface was defined as 
a function of Z and energy. The program found the closest point on 
the surface to the data point) . A two dimensional histogram of dis- 
tance from the nearest line vs. the nearest line is shown in Figure 4. 
The points in the plot on the right show the projection of this data 
on the vertical axis. Most of the data occurs between one and two 
standard deviations from the nearest line. Data close to a line occurs 
predominantly near integral lines as should be expected. Background 
events produce a slowly varing continuum, much of it in the tail of 
the curve. 

Figure 4 can be used to find the flux of individual nuclei on the 
experiment and the background. The projection of abscissa channels 
93 and 94 (in the valley between oxygen and neon) onto the ordinate 
axis is normalized upwards and is plotted as the heavy line in the plot 
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on the right. The normalization was chosen to make the areas under the 
tails of the curves equal. The similarity in shape of the tails pro- 
duced by all data and by channels 93 and 94 is excellent. In fact the 
tail has roughly the same shape at all charges. 

The peak should have the shape of a Maxwell-Boltzmann curve 17 . 

The light line in the plot at the right of Figure 4 is the sum of the 
heavy background line and a Maxwell-Boltzmann curve. The agreement 
with the data is seen to be good. 

In order to find the number of events for a given particle we inte- 
grate under the peak region of this two dimensional histogram. We then 
make corrections for the background events that lie within the region of 
integration and for the good events that lie outside the region of inte- 
gration. The integrals over the ordinate in Figure 4 for total signal and 
the background are shown .in Figure 5. The region of integration is from 
0 to 2 . 3 ct . 

We find the background under the peak at each charge by multiplying 
the tail at each charge by the peak background to tail background ratio. 
The region determined to be more than 99% background events lies above 
3.8b and our plots terminated at 9.2 a . The region which contained 
about 85% of the signal extends to 2.3a. The heavy line in Figure 4 
is certainly an upper limit to the background under the peak as some 
oxygen and neon spill into the valley and no allowance is made for 
the presence of flourine. Since a lower limit is practically zero, 
we choose the background from 0 to 2.3 a to be half of the heavy line 
and we allow the systematic error in this background to be equal to 
this background. The ratio of the background from Oq to 2.3a (peak 
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region) to the background from 3.8 cr to 9.2 <j (tail region) is found 
to be 0.286. Thus, multiplying the tail area for each charge (absissa 
bin) by 0.286 we obtain the background in the peak region given by the 
heavy line histogram in Figure 5. The shading denotes the systematic 
error. The difference between the total histogram and the background 
histogram is due to single particle events. The fraction of events in the 
tail of a Maxwell-Boltzmann curve beyond 2.3 a is 15% and is now added 
in. 


We sum counts in appropriate absissa bins for each element and 
obtain total counts and statistical and systematic errors. Next, a cor- 
rection must be applied to account for the inefficiency of the spark 
chamber and for events lost because of scattered particles that produced 
a second track in the spark chamber. 

The spark chamber efficiency was determined by finding the numbers 
of events that caused exactly 8 decks, 7 decks, and 6 decks of the 8 
deck spark chamber to fire. Assuming the average probability of a deck 


not firing is Q, the probability of having an event fire N decks is 

% l ? - N . n 

PM" (?-mV! Ni Q 


Curves of this function are shown in Figure 6. The ratios P(8)/P(7) 
and P(7)/P(6) indicate that the relative numbers of events that caused 6, 7 
and 8 decks to discharge in the particle track obey this probability formula 
reasonably well. The probability, Q, was determined for each element 
and for subgroups of events lying at various distances from the nearest 
charge line described above. The efficiency of the spark chamber is 
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different for the x and y-yiews, reflecting known readout inefficien- 
cies in the y view. A plot of the efficiencies is given in Figure 7. 

This correction is less than 10% for elements B and above in this 
flight. However, in our previous flight the chamber had only four decks, 
and with so little redundant information on the tracks Q could not be 
determined, and the separation of the primary track from its delta rays 
was more unreliable. While every attempt was made to take account of 
these factors in our previous analysis, a Q of 0.02 for a four deck 
chamber could have resulted in a 20% error in flux since we required a 
spark to be present on all four decks in both x and y. Further, this 
effect could have been Z dependant. Unfortunately, there is no way of 
determining the size or presence of this effect in our earlier data, so 
these fluxes have to stand as published 3 with the additional caveat that 
the fluxes of the elements 3 £ Z <; 10 may be systematically low due to 
small inefficiencies in the spark chamber. 

BACKGROUND 
■ Spark Chamber 

It is well known by cosmic ray experimenters that the background 
intensity falls off rapidly as Z increases away from the very abundant 
protons and alphas. This is reflected in our data also. In Figure 8, 
the probability of a second track being registered in the spark chambers 
is shown to decrease as the charge increases. This probability goes 
from a few percent at Li to below 1% at C and leads us to believe that 
trajectory confusion due to backscattering events is not a significant 
source of error. This is further verified by Table 1 where it is seen 
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that well identified nuclei are correlated with good spark chamber 
tracks. The background of well identified events with no spark chamber 
tracks is probably due to air showers, and falls off rapidly with 
increasing pulse height. This background makes the identification of 
Li very difficult. 

From this table we see that any procedure which places a requirement 
that a good track be present greatly enhances the signal to noise ratio 
and does not significantly bias against good events. 

The atmosphere and the material in the experimental apparatus 
cause cosmic rays to spallate and lose energy . Interactions that caused 
two or more spatially resolvable particles to pass through the spark 
chamber were looked for in the data. If two or more tracks diverged from 
a point in the x-view and the vertical height of this point from the 
center of the spark chamber equalled the vertical height of a similar 
point in the y-view, an interaction was considered to have occurred. 

One such event is shown in Figure 9. We plot the interaction vertices 
of all such events in both the x and y-views and obtain the scatter plot 
in Figure 10. The locations of the detectors are indicated and the inter 
actions can be seen to originate in materially dense regions. The number 
of these events that caused at least 7 or 8 decks to discharge in each 
of two tracks is 633 or 4% of the single ion events. When the criterion 
of a good fit was applied to these events they became a negligible 
fraction. The probability of the spark chamber registering such an 
event is approximately the square of the probability given in Figure 7 
or about 0.8 above B. In the analysis of the data these interaction 
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events were not explicitly considered. Most of them interacted in the 
charge module and were accounted for by (1) the background under the 
Maxwell-Boltzmann curve and (2) the final interaction cross section 
correction. Perhaps as many as 100 events possessed backs cat tered 
interactions from the top of the spectrometer. This is less than 1% 
of the total events and seems to be an unlikely source for any large 
systematic error. 

The correction for the efficiency of the spark chamber to have 
more than one track in a view was made assuming the additional tracks 
were delta rays and hence these events should be included. If some of 
the multiple track events were actually spallation events, the correc- 
tion should not be as large as given in Figure 7. The probability of 
spallation increases with charge. Figure 7 shows the correction to be 
0.5% for high charges and so we have overestimated our fluxes at most 
by this amount. 

Maxwellian Distribution 

In order to more fully understand the background in the tail of 
the curve in Figure 4, we looked at the pulse heights of all detectors 
for 68 carbon events satisfying (1) 33 £ S 1 ^ 43 and (2) distance 
from nearest charge line ^5. We made 5 catagories of events and found 
the following distribution: 

1. 46 (68%) events characterized by 

a. one detector having an amomalously low or high signal and 

b. the following detectors having anomalously low signals. 

These events are probably interactions. 
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2. 5(7%) in which the pulse height was anomalously high, and all 
other pulse heights were normal. These could be due to a delta 

ray going through a PMT tube face. (This effect is unimportant in C^.) 

3. 8(12%) in which one scintillator is too high. These fluctuations are 

at least 5<j and are unlikely to be caused by Landau fluctuations. 

They could be due to nuclear excitations of detector atoms. 

4. 5(7%) spectrometer signal anomalously large due to the non- 

Gaussian distribution of signals in the thin (0.75 mean free path) 
spectrometer. 

5. 4(6%) showers and other explanations. 

Thus a detailed look at events in the tail of the goodness of fit curve 
shows that most involve interactions. 

In order to account for the possibility that events in categories 2-4 
are legitimate, we have included a correction for good events under the 
background curve from 2.3cr to 9.%y. We have linearly added 50% of these 
added events to the systematic error. This correction accounts for a varying 
amount of the flux. For example the correction is 5% for carbon and 25% for 
a rare element like Na. This represents the least satisfying part of our 
data analysis. 

Category 2 could have been eliminated if we had divided the Cerenkov 
phototubes for each detector into two sets and checked for consistency 
between them. 

A correction could have been made for Category 3 by knowing the excitation 
cross sections, but this category of error is probably common to all compo- 
sition measurements made with plastic scintillators. 

Category 4 could have been eliminated by the use of a thicker spectrometer. 
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Corrections 

Corrections for spallation in the atmosphere are made for the aver- 
age 3.91 g/cm 8 average depth during the data taking period and for the 
7.25 g/cm 2 of experiment thickness through the charge determining module. 
The mean free paths for the various elements were obtained from Webber 
et al 14 . No errors from uncertainties in the mean free paths are included 
in our analysis. 

A bias was inadvertently introduced in one of our spark chamber 
track selection criteria that primarily affects the lowest elements. 

Since the errors in these fluxes are already very large, we have simply 
corrected Li upward by 10%, Be by 6%, B by 4%, C by 1% and with no 
corrections for heavier elements. We have linearly increased the res- 
pective systematic errors by 3%, 2%, 1% and 0% of the respective abun- 
dances. We felt it unjustified to redo our analysis because the results 
would not be appreciably changed. 

Results 

The geometrical factor of the 50 cm square by 80.5 cm high experi- 
ment is 776 ctt^ star. Live time accounted for 55.67% of the 19 hours 
and 19 minutes of the data taking period yielding a live exposure time 
of 38713 seconds. The raw data and all the corrections applied are shown 
in Table 2. The symbols used in Table 2 are defined in Figure 11. 

The resulting fluxes of the elements and the statistical and system- 
atic uncertainties are given in Table 3. The results of a complete 
analysis without including the spectrometer data is also given. These 
results are plotted in Figure 12. Various ratios are given in Table 4. 

The systematic errors play a smaller role in the uncertainty of a ratio 
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than in the flux because the errors in the numerator and denominator both 
vary in the same way. The systematic error in the ratio is obtained by 
assuming the background is either underestimated or overestimated in both 
elements. The error of a ratio in Table 4 is the square root of the sum 
of the square of the fluctuation due to systematic errors and the square 
of the statistical error. 

It can be seen in Figure 12 that the background in the region of 
the light nuclei is substantially reduced by including the spectrometer 
in the analysis. Also note the effect of the spectrometer in lifting 
the degeneracy in response between high energy particles of a low charge 
and lower energy particles of the next higher charge. One might expect 
the sysyematic errors in the background correction to produce uniformly 
high or low abundances for a given analysis . That they do not suggest 
that the background distributions under the peak and/or the contributions 
of good events in the tail are not known. They are probably coupled 
to the analysis method, differently for different charges. It further 
reflects that a careful analysis of systematic errors is essential 
to making interexperiment comparisons of absolute intensities. 

COMPARISON WITH PREVIOUS DATA 

In Table 4 our flux ratios are presented and compared with previously 
published values. 15 * 4 The relative abundances are in generally good 
agreement, but there are still some differences that are outside the quoted 
systematic errors, even if we assume other experiments have systematic 

errors as large as ours . 

Previous carbon and oxygen intensity measurements and our data point 
are shown in Figure 13. There is agreement with Cal Tech 19 and disagree- 
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ment with the Universities of New Hampshire and Chicago . In order to try 
to understand the disagreement, we analyzed the first ~ 10,000 events in our 
data with only one selection criterion: that (S^ - S^/S^ + S 2^ ^ This 

criterion is very roughly the criterion used in the experiments of the New 
Hampshire and Chicago groups which are similar to ours . Neither of these 
experiments contans a spark chamber, and so if no track could be found, the 
particle was assumed to be vertically incident. The resulting histogram vs. 
distance from the nearest charge line is given in Figure 14. A background 
cruve obtained from the flourine valley and a curve of this background plus a 
Maxwell-Boltzmann distribution are also shown. We find that ~ 1/3 of the 
events that satisfy 2(S^ - S 2 )/S^ + S^) are background events according to 
our analysis. 

This may be in part responsible for the differences between intensities. 

The analysis conducted by the Cal Tech group gives an intensity that agrees 
with ours. Their telescope is made of solid state detectors and the systematic 
effects should be the least. Their spectrum, however, is flatter than that 
of the New Hampshire group at low energy. 

In summary, we believe that the discrepancy between the fluxes of 

2 X2 4 

Balasubrahmanyan and Ormes and those of Webber et al. , Julliuson and Smith 

20 

et al. is partially due to spark chamber inefficiencies in the former experi- 
ment and partially to background rejection criteria in the latter experiments. 
It Is important to resolve this discrepancy, because it affects the interpre- 
tation of the C+O/Fe ratio which is dependent upon energy. We have found 
nothing in this analysis which would lead us to doubt our earlier result 
that the Fe spectrum is significantly flatter than the C+0 spectrum. What 
is at issue is the relative intensities. If the C+O/Fe ratio at earth can 
be shown to go below the calculated source ratio at energies above 10 GeV/arau, 
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then there is no extreme energy dependent propagation model which can 
explain it, and one must turn to energy dependent models at the source 
or different sources of different species of cosmic rays. 

We hope to resolve this deiscrepancy by flying a spectrometer 
sufficiently deep to measure energy along with a gas Cerenkov 
detector to calibrate the spectrometer for Z > 1 nuclei at 15 GeV/amu. 
The data will then be analyzed using the techniques described here. 
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FIGURE CAPTIONS 


Figure 1 
Figure 2 

Figure 3 

Figure 4 

Figure 5 


Figure 6 

Figure 7 
Figure 8 

Figure 9 


Experiment configuration. 

Scatter plot of average scintillator pulse height vs. average 
Cerenkov pulse height. Solid lines are drawn by eye. 

Scatter plot of summed spectrometer pulse heights vs. average 
Cerenkov pulse height for oxygen. Solid line is drawn by eye 

Scatter plot of distance from nearest charge line vs. nearest 
charge. Plot at right is projection of data on vertical axis 
Heavy line is background. Light line is background plus a 
Maxwell-Boltzmann curve normalized to peak data point. 

Signal and background. The heavy line is background; the 
shading shows the range of the systematic error. The signal 
as shown by the light line is coincident with the background 
between low elements but is higher than background between 
higher element peaks due to the wider tails. 

Probability of a track containing N decks with discharges in 
the track in an 8 deck spark chamber vs. the probability of 
a deck not discharging in the track. 

Probability of a particle track discharging 7 or 8 decks. 

Probability of a particle track that discharged 7 or 8 decks 
to be accompanied by another 7 or 8 discharged deck track. 

X and Y views of spark chamber that recorded event with two 
tracks. Each dot represents two wires in the wire grid spark 
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chamber. "L" denotes the left hand wire discharged; "J" the 
right hand wire; and "U" both wires. 

Figure 10 Coordinates of interaction vertices in X and Y-views. Inter- 
actions occur where massive parts of the experiment lie. 

Figure 11 Definition of regions under the Maxwellian curve. 

Figure 12 Abundances with and without using the spectrometer data in 
the analysis. 

Figure 13 Previously measured carbon and oxygen abundances and our data 
point. 


Figure 14 Number of events vs. distance from nearest charge line for 
events that satisfy (Sj - S^)/(S t + ) £ 0.3. The expected 

Maxwellian curve and the background, given by scaling up the 
flourine valley, are also given. 



TABLE 1 

EFFICIENCY TABLE 


S 1+S 2 

Distance 

from 

line in 

units of cr 



2 








0-2 

2-4 

4-8 

None 




337 

605 

1199 

8 

0 

Number 

0-18 

213 

168 

180 

57 

1 

of 


14 

14 

17 

1 

^2 

tracks 


Distance 

from 

line 





0-2 

2-4 

4-8 

None 




78 

156 

1074 

8 

0 

Number 

18-36 

271 

178 

288 

0 

1 

of 


2 

10 

22 

0 

>2 

tracks 


Distance 

from 

line 





0-2 

2-4 

4-8 

None 




32 

50 

581 

5 

0 

Number 

36-54 

431 

278 

271 

0 

1 

of 


3 

1 

20 

0 

Ss2 

tracks 


Distance 

from 

line 





0-2 

2-4 

4-8 

None 




22 

28 

344 

9 

0 

Number 

54-72 

293 

173 

156 

13 

1 

of 


2 

1 

13 

0 

>,2 

tracks 


Distance 

from 

line 





0-2 

2-4 

4-8 

None 




6 

16 

504 

34 

0 

Number 

72-190 

243 

168 

282 

35 

1 

of 


2 

0 

15 

1 

£2 

tracks 


The total data consists of 168527 events. The first 10000 events were 
analysed and 714 tracks were found that passed outside C-^ . 8922 of the 
remaining events have (S^+S 2 ) /2 < 190 and are included in this table. 

Each track must discharge 7 or 8 decks. The column labeled "none" contains 
events that lie closest to the helium or iron lines or that have multiple 

minima. Multiple minima can occur if an event lies very far from the res- 
ponse surface. Note that no selection on the consistency of Sj. and % has 


been made. 



TABLE 2 



bins 

No-2.29 a 

Na 

p sig 

Pstat 

p sys 

Q sig 

°stat 

ttsys 

*1 

e 2 

ix. 

g.f. 

Flux 

AF s tat 

£ 

LL. 

c 

Li 

24-33 

1016 

374 

1224 

38 

192 

1346 

42 

232 

.74 

.948 

.796 

3004 

.798 

.025 

.138 

Be 

37-43 

574 

236 

693 

26 

121 

735 

30 

135 

.86 

.961 

.775 

3004 

.379 

.015 

.070 

B 

47-54 

1526 

528 

1B36 

46 

272 

1909 

48 

291 

.91 

.981 

.755 

3004 

.943 

.023 

.144 

C 

55-65 

5018 

775 

5977 

84 

480 

6037 

85 

400 

.925 

.991 

.687 

3004 

3.191 

.045 

.211 

N 

66-74 

1438 

710 

1744 

45 

366 

1744 

45 

366 

.93 

.995 

.686 

3004 

.915 

.024 

.192 

0 

d 

75-69 

4401 

693 

5243 

78 

358 

5243 

78 

358 

.935 

.995 

.863 

3004 

2.831 

.042 

.193 

r 

Ne 

94-106 

742 

259 

893 

32 

134 

893 

32 

134 

.94 

.995 

.653 

3004 

.487 

.017 

.073 

Ne 

107-113 

202 

110 

245 

17 

56 

245 

17 

56 

.94 

.995 

.664 

3004 

.131 

.009 

.030 

Mg 

114-127 

785 

134 

936 

33 

69 

936 

33 

69 

.94 

.995 

.821 

3004 

.537 

.019 

.040 

A 

128-132 

117 

66 

142 

13 

34 

142 

13 

34 

.94 

.995 

.631 

3004 

.080 

.007 

.019 

Si 

133-147 

674 

98 

802 

31 

51 

802 

31 

51 

.94 

.995 

.590 

3004 

.484 

.019 

.031 


bitu denotes bins along abscissa in Figure 5 
Nq_2.29<f» N|j, N c are defined in Figure 11 
*} = N b /N fl = 0.273 

f = N c /N, = OSS 

G = Number of good events in tail = 0.25 N fl (1 + f) = 0.3675 N g 

Ojj * Area under Maxwellian eunre/Area under Maxwellian curve from 0 to 2.29 cr * 1 .1 81 

P si0 = Number of good events = <N 0 _ 2 2 g a - fjNg) g„ + G = 1,181 N fl _ 2 2 g c + 0.065 N a 

P sta , - Statistical error = J 1.1 81 2 N 0 _ 2 29{J + 0.065 2 « 8 * J7.395 N 0 -2.29o + 00042 N a 

Py S = Systematic error = r?g K N g + 0.5 G = 0.516 N a 

Qjjg hra correction for analysis error. = P^x (1.10 for Li, 1.06 for Be, 1.04 for B, 1.01 for C, 1.00 for others) 
Q Jtat has correction for analysis error. Q Jtat = P (to1 x (1.10 for Li, 1.06 for Be, 1.04 for B, 1.01 for C, 1.00 for others) 
°sys hes correction for analysis error. Q^, * + Q^. x (0.03 for Li, 0.02 for Be, 0.01 for B, 0.0 for others) 

= Efficiency for one or more tracks in spark chamber 

= 1 minus efficiency for two or more tracks in spark chamber 

i.c. = interaction correction 

B-f. - Geometrical factor 

Flux = Ojig/ e l e 2 (ix.) (g.f.) 

4F stat s <W e 1 e 2 < ic ^8 f ) 

^ F sys ~ ®sys ^ ^ 1 ^ 2 ®’*.) (g.f.) 





TABIE 3 


Element 


With Spectrometer 



Without Spectrometer 



Flux 

Errors 


Flux 

Errors 




Statistical Systematic 


Statistical Systematic 

Li 

0 82 

.03 

.15 

00 

. 

1— 1 

.03 

.10 

Be 

.39 

.02 

.07 

.50 

.02 

.10 

B 

.96 

.02 

.15 

1.00 

CM 

o 

0 

.21 

e 

3.22 

.05 

.23 

3.35 

.04 

.23 

N 

.94 

.02 

.21 

.84 

.02 

.20 

0 

2.86 

.04 

.21 

2.71 

.04 

.26 

F 







Ne 

.50 

.02 

.08 

.43 

.02 

.09 

Na 

.14 

.01 

.03 




Mg 

.54 

.02 

.04 

.52 

.02 

.06 

Al 

.082 

.007 

.02 




Si 

.49 

.02 

.03 

.53 

.02 

.04 


n 

Units are number/m sec ster 



TABLE 4 


Ratio 


of 

With 


Fluxes 

Spectrometer 

Be/B 

0.40 

+ 

.02 

Li/O 

0.28 

± 

.03 

Be/O 

0.135 

+ 

.02 

B/0 

0.34 

+ 

.03 

c/o 

1.13 

± 0 

.02 

N/O 

0.325 

+ 

.05 

Ne/O 

0.17 

+ 

.02 

Na/O 

0.047 

+ 

.009 

Mg/O 

0 . 190 

+ 

.007 

Al/O 

00 

CM 

O 

o 

± 

.006 

Si/O 

r-t 
t". 
T— l 
0 
o 

+ 

.007 

Li+Be+B/C-W+O 

0.31 

+ 0 

.03 


Webbe 


et al 15 
> lGeV/amu 

Julius son^ 
@6 .0gm/cm 
1 GeV/amu 

.293 

± *010 

0.38 

.192 

± - 004 

0.27 

.101 

+ .003 

0.15 

.345 

± *005 

0.39 

1.13 

+ .02 

1.17 

0.286 

+ .005 

0.33 

0,174 

+ .006 

0.16 

0.030 

+ .002 

0.04 

0.204 

+ .005 

0.20 

0,030 

+ .002 

0.04 

0.137 

+ .004 

0.15 

0.265 

+ .01 

0.324 
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